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ABSTRACT 

We study the redshift evolution of the luminosity function (LF) and redshift selec¬ 
tion effect of long gamma-ray bursts (LGRBs). The method is to fit the observed peak 
flux and redshift distributions, simultaneously. To account for the complex triggering algo¬ 
rithm of Swift, we use a flux triggering efficiency function. We find evidence supporting 
an evolving LF, where the break luminosity scales as L* oc (1 h- zY, with r = 3.5 ^q 2 and 
T = for two kind of LGRB rate models. The corresponding local GRB rates are 

R(0) = 0.86;^Qggyr“'Gpc“^ and R(0) = 0.54;^°yr * respectively. Furthermore, by 

comparing the redshift distribution between the observed one and our mocked one, we find 
that the redshift detection efficiency of the flux triggered GRBs decreases with redshift. Es¬ 
pecially, a great number of GRBs miss their redshifts in the redshift range of 1 < z < 2.5, 
where “redshift desert” effect may be dominated. More interestingly, our results show that the 
“redshift desert” effect is mainly introduced by the dimmer GRBs, e.g., P < 10“^ erg/ s/ cm^, 
but has no effect on the brighter GRBs. 
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1 INTRODUCTION 

Gamma-ray bursts (GRBs) are one of the most luminous and dis¬ 
tant transients in the universe (Greiner et al. 2009; Cucchiara et 
al. 2011). Observationally, GRBs are usually categorized into two 
groups: spectrally soft long GRBs with Tgo > 2s are expected 
to result from the collapse of short lived massive stars, as some 
long GRBs are evidenced to connect with Type Ic supemovae (e.g. 
Galama et al. 1998; Bloom et al. 2002; Stanek et al. 2003; Thomsen 
et al. 2004; Campana et al. 2006; Berger et al. 2011; Melandri et al. 
2012). Meanwhile, spectrally hard short GRBs with Tgo < 2 s are 
believed to originate from the merger of compact stars (e.g. Eichler 
et al. 1989; Nakar 2007; Fong et al. 2013). In this paper, GRBs are 
actually referred to long GRBs, unless otherwise specified. 

The hosts of GRBs in the star formation regions suggest that 
GRBs could be used as the tracer of star formation rate (SFR) 
(Totani 1997; Paczyhski 1998; Wijers et al. 1998). The detection 
of GRB 980425 associated with a supernova also strengthened the 
expectation (Galama et al. 1998). Therefore, GRBs provide a new 
opportunity for the measuring of star formation history (for a re¬ 
cent review, see Wang, Dai & Liang 2015), especially at high red¬ 
shift, where direct measurement is difficult. Thanks to the launch 
of Swift satellite, which provides a large number of GRBs with 
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measured redshifts. This large GRB sample makes it possible to 
give a more tight constraint on the relation between GRB rate and 
SFR, and rules out the models that GRBs unbiased trace the cos¬ 
mic star formation history (Chary et al. 2007). Recent results show 
that GRBs do not trace the star formation history exactly but with 
an additional evolution, e.g., (1 -l- zY (Le & Dermer 2007; Kistler 
et al. 2008; Wang & Dai 2009; Wanderman & Piran 2010; Cao et 
al. 2011; Robertson & Ellis 2012; Wang 2013). However, the value 
of 6 varies large from 6 = 0.5 to 5 = 2 which strongly depends on 
the sample selection, and the inferred SFR at high redshift seems 
too high comparing with the observation from the galaxy surveys 
(Kistler et al. 2008, 2009; Robertson & Ellis 2012; Wang 2013). 
Furthermore, Yu et al. (2015) even found an excess of GRB rate 
at low redshift of z < 1.0 (see also Petrosian et al. 2015). Many 
theoretical models have been proposed to explain the discrepancy 
between the GRB rate and SFR, e.g., the evolution of the cosmic 
metallicity (Danger & Norman 2006; Li 2008; Elliott et al. 2012), 
the evolution of the initial mass function (Wang & Dai 2011), or 
the additional cosmic string explosions (Cheng et al. 2010). How¬ 
ever, most previous works are based on the assumption that the lu¬ 
minosity function (LF) is a constant form and independent of red¬ 
shift. Alternatively, some works suggest that the LF should evolve 
with redshift by using the GRB sample with pseudo redshifts in¬ 
ferred from luminosity relation (e.g., luminosity-variability rela¬ 
tion, luminosity-peak energy relation; Lloyd-Ronning et al. 2002; 
Firmani et al 2004; Yonetoku et al. 2004; Matsubayashi et al. 2005; 
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Kocevski & Liang 2006; Tan et al. 2013). The evolution of the LF 
scales as L oc (1 + zT, with t varies from 1 to 3. In this paper, we 
will use GRBs with measured redshifts, but not the pseudo ones, to 
find out whether LF evolves with redshift or not. 

In this work, we specially pay attention to two important ef¬ 
fects. One is the redshift selection effect. Previous works use the 
GRB redshift distribution to study the LF or GRB rate do not con¬ 
sider the missing redshift problem (e.g. Kocevski & Liang 2006; 
Salvaterra et al 2012). However, this selection effect is very impor¬ 
tant for these studies, as only ~ 30% of Swift GRBs have measured 
redshifts. Especially, if the selection effect evolves with redshift, 
the results even could be wrong. A detailed study about the red¬ 
shift selection effect is presented in Coward et al. (2013) (see also 
Fiore et al. 2007). The other effect is the Swift triggering thresh¬ 
old problem. Before Swift, a single detection threshold based on an 
increased photon count rate above the background is a reasonable 
approximation. However, Swift has a much more complex trigger¬ 
ing algorithm in order to maximize the GRB detection (Band 2006), 
e.g., the Burst Alert Telescope (BAT) on board Swift has over 500 
rate trigger criteria. Therefore, a single detection threshold approx¬ 
imation may not be correct for Swift. To solve this problem, we use 
a more complex triggering algorithm which is based on the most 
recent result of Lien et al. (2014). 

First, we study the redshift evolution of the LF by fitting the 
number distributions of peak flux and redshift. Considering the dis¬ 
crepancy between the SFR and GRB rate, we introduce two kind 
of GRB rate models. While fitting the redshift distribution, only 
redshifts measured by absorption spectroscopy are considered (to 
reduce the redshift selection effect; see the following section for 
details). For the peak flux data, we take it from the Butler’s on¬ 
line catalog Q Then, we analyze the redshift selection effect by 
comparing the number distribution between the observed redshift 
sample (include all GRBs with measured redshifts) and the model 
predicted one. We find that the selection effect is not only redshift 
dependent but also flux dependent. 

This paper is organized as follows. In section 2, we describe 
the data extraction methodology. In section 3, models applied to 
constrain the LF and redshift selection effect are presented. Results 
are presented in section 4. The summary and discussion will be 
given in section 5. 


2 THE DATA 

Including long and short GRBs, Swift has detected more than 900 
GRBs (until GRB 141026A) in the past 10 years. However, recent 
studies suggest that the classical classification of short and long 
GRBs with Tgo = 2 s is not appropriate for Swift. Bromberg et al. 
(2013) argued that a more suitable selection for long GRBs from 
Swift satellite should be defined by Tqo > 0.8 s, which is based 
on a physically motivated model. Here we believe that the clas¬ 
sification with Tgo = 2 s is more strict in our study. As shown 
in Figure 3 of Bromberg et al. (2013), the probability of a GRB 
with Tgo > 2 s to be a collapsar is more than 80%, which is higher 
than 50% for Tqo > 0.8 s. Additionally, three sub-luminous GRBs 
(GRB 060218, 060505, 100316D) with luminosity L < 10'**’ erg s^* 
have been excluded from our sample (Soderberg et al. 2006; Liang 
et al. 2007; Virgilii et al. 2008; Howell et al. 2011; Howell & 
Coward 2013), three ultra-long GRBs (GRB 101225A, 111209A, 
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Figure 1. Top: The redshift distribution of the full GRB sample including 
244 GRBs. Bottom: The reduced redshift distribution of 182 GRBs with 
redshifts obtained either from only afterglows, or from both afterglows and 
host galaxies. 


121027A) have been suggested to be another population(Gendre et 
al. 2013; Stratta et al. 2013; Levan et al. 2014). We also exclude 22 
possible short GRBs (GRB 050416A, 050724, 050911, 051016B, 
051114, 051227, 060505, 060614, 061006, 061210, 070714B, 
071227, 080123, 080503, 080520, 09053IB, 090607, 090715B, 
090916, 100213A, 100216A, 100816A; Dietz 2011; Kopac et al. 
2012; Fong et al. 2013; Tsutsui et al. 2013; Howell & Coward 2013; 
Berger 2014; Howell et al. 2014). 

In this paper, the redshift data was taken from the Swift archive 
0. After removing the uncertain redshifts, we obtained 282 GRBs 
(including long and short GRBs) with measured redshifts, which is 
around 30% of the BAT-triggered GRBs. For long GRBs, their red¬ 
shifts are always measured from afterglows and/or the host galax¬ 
ies observed by the ground-based telescopes, including the Very 
Large Telescope (VLT), Gemini-S-N, Keck and Lick. Finally, we 
obtained 244 long GRBs with confirmed redshifts (the full redshift 
sample). However, this full redshift sample is greatly affected by 
the observation biases. On one hand, more low redshift GRBs could 
be observed, as the brightest GRBs are predominantly nearby. On 
the other hand, some redshift dependent selection effects even 
could change the shape of the redshift distribution. As shown in 
the upper panel of Figure 1, a large number of GRBs missed their 
redshifts in the range of 1.5 < z < 2, which may caused by the 
“redshift desert” effect. Therefore, to derive a less biased redshift 
sample is an immediate necessity of this work. Following Fynbo et 
al. (2009), we reconstruct our redshift sample by considering red¬ 
shifts measured either from afterglows, or from both afterglows and 
host galaxies. We exclude the bursts with redshifts measured from 
the emission spectra of the host galaxies, as these bursts always 
cover a low redshift range (e.g., 0.3 < z < 2.8; Wanderman & Piran 
2010, Howell et al. 2014). We also exclude GRBs with photomet¬ 
ric redshifts because of the large uncertainties. After doing this, we 
obtained a less biased redshift sample with 182 GRBs (the reduced 
redshift sample). The distribution of the reduced redshift sample is 
shown in the bottom panel of Figure 1. It’s obvious that the com¬ 
pleteness of the reduced redshift sample is better than that of the 
full redshift sample, as shown in the top and bottom panels of fig¬ 
ure 1 (e.g., in the redshift range of 1.5 < z < 2). However, one may 
argue that the missing redshifts of GRBs in the range of 1.5 < z < 2 
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Figure 2. Top: The triggering efficiency of Swift as a function of peak flux 
(blue solid line), and the single detection threshold of 2.6 ph s“* cm“^ (gray 
dashed line) and 0.4 ph s“* cm“^ (gray dotted line). The middle and hottom 
figures show the effects of the triggering algorithm on the number distribu¬ 
tions of peak flux and redshift (for RGRBl for example), respectively. 


may be caused by the statistical fluctuations. For this point, we will 
discuss it in section 4.1. 

In this work, we use the GRB peak flux data in the energy 
range of 15-150 keV of the BAT. The data was taken from the But¬ 
ler’s online catalogue, which is the extension of the work presented 
in Butler et al. (2007, 2010). As the peak flux data can be measured 
directly from BAT by assuming the fewest burst characters, we sug¬ 
gest that the peak flux distribution is the least uncertain property. 
Finally, we obtained a total number of 681 bursts with peak fluxes. 


3 THE MODEL 

The first step of our work is to fit the GRB number distributions 
of peak flux and redshift, which are both related to the GRB rate 
and LF. The expected number of GRBs with observed peak fluxes 
between and P 2 that triggered the BAT can be expressed by 

AO r°° F- ■ dV(z) 

N{PuP 2)=—T\ rj{P)^,Am{z)dL-^, (1) 

47r Jo Jij l+z 

where ~ 0.1 is the field view of the BAT, T ~ 10 yrs is 

the observational period, dV(z) is the comoving volume and 1/(1 -f 
z) accounts for the time dilation, rj(P) is the triggering function, 
d?.p{L) is the LF with = lO'*^ erg s“* and = 10^^ erg s“' 
for normalization, and P(z) is the GRB rate. The expected number 
of GRBs within redshift range of < z < Z 2 is given by 

Afi r- r'"”“ . dV(z) 

N(Zi,Z2)=—T r,(P)<t>,AQRiz)dL-^. (2) 

47r J-, 1 + z 

For the flat ACDM cosmology, we employ the cosmological pa¬ 
rameters from WMAP nine-year results with Q„, = 0.28, Aa = 0.72 
and Ho - 70 km s“* Mpc^'. 

The triggering algorithm of Swift is complicated. Lien et al. 
(2014) simulated 50,000 GRBs to mock the triggering efficiency. 


where GRB light curves, incident angles, BAT’s active detector 
number, and the image trigger algorithm are considered. Although 
a specific set of parameters are included, we consider that the trig¬ 
gering algorithm is better than other techniques. By comparing the 
simulated peak flux number distribution with the real triggered one, 
Howell et al. (2014) derived the flux triggering efficiency in a func¬ 
tional form of 


u(R) = 


a(b + cPIPq) 
(1 -h P/dPo) ’ 


(3) 


with P > 5.5 X 10“^ erg s“' cm“^, and below this range the function 
equals to zero. The function parameters are as follows: a=0.47, b=- 
0.05, c=L46, d=L45 and Pq = 1.6 x 10“^ erg s“* cm“^ (Howell 
et al. 2014). We could reproduce the triggered GRB population in 
Lien et al (2014) by using this triggering function. It seems that the 
sample completeness is much better than that of the step function 
approximation. In figure 2, we compare t]{P) with two step func¬ 
tion approximations, e.g., a value of 0.4 ph s“* cm“^ and a value of 
2.6 ph s“* cm“^, which have been used by Salvaterra et al. (2012). 
From the middle panel of the figure, we see that the threshold of 
0.4 ph s ' cm“^ predicts more dim GRBs, while the threshold of 
2.6 ph s ' cm“^ predicts more bright ones. It is apparent that the 
adoption of such approximations could be problematic while esti¬ 
mating the relative contributions of bright and dim GRBs. Further¬ 
more, it seems that the adoption of different triggering algorithms 
do not change the shape of the redshift distribution significantly, but 
change the GRB number greatly as shown in the bottom panel of 
the figure. So it is obvious that the threshold of 0.4 ph s“* cm“^ pre¬ 
dicts a lower local GRB rate, while the threshold of 2.6 ph s“* cm“^ 
predicts a higher one. 

The isotropic-equivalent peak luminosity in the source frame 
can be calculated by 




L = ATTdi^izfP^, 
b 


(4) 


where P is the peak flux (in units of erg s“' cm“^) in the observed 
energy band of 15 - 150 keV, dz{z) = (1 +z)rfc(z) is the luminosity 
distance. The term k(z)lb is to convert the observed energy band to 
the rest frame band of 1 - 10'* keV. Here b accounts for the bolo- 
metric energy fraction that seen in the detector band (Wanderman 
& Piran 2010, Howell et al. 2011), and k(z) is the cosmological 
correction. We express them in the form of 


pl50 

b = j ES(E)dEI I ES(E)dE 


1 


and 


pl50 p 

k(z)= ES{E)dEI 

Jl5 Jlf 


150(l+z) 


ES(E)dE, 


I5(l+z) 


(5) 


(6) 


respectively. Here S (E) is the rest frame photon spectrum, which 
can be well expressed by the empirical Band function (Band et al. 
1993, 2003). The high and low energy spectral indices are given by 
-2.25 and -1, respectively. The spectral peak energy in the source 
frame can be derived by the Yonetoku relation (Yonetoku et al. 
2004) 


Edl+z) = 


12.34 X 10^’j 


(7) 


which is a much tighter and reliable relation. 

We assume a broken power law form for the LF, which is sug¬ 
gested to be better than that of a single power law(e.g., Cao et al. 
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Figure 3. Left: The number of peak flux distribution (top; PFD) and redshift distribution (bottom; RD) of RGRBl. The solid blue lines are our best fitting 
results with an evolving LF (ELFl). The dashed gray lines are the fitting results with an non-evolving LF (NLFl). Right: The same as the left figures but for 
RGRB2. 


2011; Tanetal. 2013) 


(2010), we describe it in a functional form of 


[ 

1 / f \ 

L < Lt(z}, 

( 8 ) 

Mz)=Ro\ 

1 (ilbi) ’ 

L>Lt(z), 


( 


(l+z)’'‘, 

(l+^)’'^ 


z < z„ 
z, < z. 


( 11 ) 


where the break luminosity is assumed to evolve with redshift in the 
form of Li(z) = C x (1 -l- z)'^. The normalization coefficient is taken 
by assuming the minimum luminosity of L,nin = 10“*^ ergs“*. Here 
a, /3, T and C are the free parameters which will be determinated in 
the following section. 

Finally, we consider two kind of GRB rate models. For the 
first one, we consider that GRB rate follows the SFR (RGRB1 for 
short). The observational GRB production rate can be connected to 
star formation rate as 

Rx(z) = McP.iz), (9) 

where ft, = 0.01 is the beaming degree of GRB outflows and 
/(, arises from the particularities of GRB progenitors (e.g., mass, 
metallicity, magnetic field, etc). We assume ftfc does not evolve 
with redshift, and the redshift evolution effect between GRB rate 
and SFR mainly originates from the LF. 

For the SFR p,(z) , we describe it as (Hopkins & Beacom 
2006) 

r (l+z)^■'*^ z<0.97, 

p.(z)«j (1 + z)”, 0.97 < z < 3.5, (10) 

[ (l-i-z)L 3.5 <z, 

with the local star formation rate p,(0) = 0.02Moyr“*Mpc“^. Here 
we assume the high-redshift SFR (z > 3.5) evolves in the form of 
a power law, which is still ambiguous. However, motivated by the 
shapes of the SFR at z < 3.5 and also implied by some preliminary 
measurements (Bouwens et al. 2008, 2011; Ellis et al. 2012) and 
GRBs as well as their host galaxies (Yiiksel et al. 2008; Kistler et 
al. 2009; Wang & Dai 2011; Elliott etal. 2012), the power low form 
of the SFR is considerable. For simplicity, we set a typical value of 
K = -3 in the following calculation (e.g., Wang 2013). 

The second GRB rate model (RGRB2 for short) is derived 
from the observed long GRBs. Following Wanderman & Piran 


with typical values of z, = 3.6, Vi = 2.1, and V 2 = -0.7 based on 
the recent study (e.g., Howell et al. 2014). This form of GRB rate 
model is widely used in literature. 


4 RESULTS 

In this section, we constrain the model parameters by jointly fit¬ 
ting the peak flux distribution (PFD) and redshift distribution (RD; 
the reduced redshift sample). As mention in the above section, two 
GRB rate models are considered. In each model, five free param¬ 
eters are to be constrained: a, r, C, and (for RGRBl) or Rq 
(for RGRB2). In fact, fc or Ro can be easily derived by equating the 
model predicted peak flux number (equation 1) to the observed one 
(number of 681), if other four parameters are constrained. Other 
three parameters of a, and C are mainly determinated by the 
PFD. The lower P-values of PFD fix the values of a and C, and 
fj is mainly determined by the higher P-values. For the value of r, 
it is strongly dependent on the redshift evolution of RD. 

Before the start of the work, we first set a primary constraint 
on the related parameters, e.g., a > 0, /? > 0, L/,(0) > lO"*® erg s“', 
r > 0, and the local GRB rate R(0) 6 [0,2] yr“* Gpc“^ for both 
GRB rate models. Then, we give an arbitrary set of values for the 
five free parameters in the allowed ranges. For each set of the pa¬ 
rameters, we calculate the;^'^ values for both PFD (;yppQ) and RD 
The totalis assumed to be the linear combination of;^fppJJ 
and Xrj}- The best fit parameters are derived by minimizing the 
global At last, we obtain the error bars for the five parameters 
by setting = x^in 5.89, which corresponds to the parameters 
within 68.3% confidence level. Table 1 gives the values of the fitted 
parameters for each model. The degrees of freedom (for the simul¬ 
taneous fit) in Table 1 is 28 minus the number of parameters. The 
information about the quality of the fit for models is also given in 
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1 q+0.05 
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d.j_0_2 

0-33ti55 

4 04+0.6 
-0.39 

o-86l:‘‘ 

17.74 

0.77 

12.31 

0.26 

5.43 

0.71 

ELF2 

Q j+0.73 

9 1 +0.13 
-0.03 

^•°-0.08 

141 

— 

0.541:^5 

23.63 

0.42 

16.43 

0.09 

7.20 

0.51 

NLFl 

0 

^•^-0.1 

1 "7+0.05 

24^8:?^ 

^•^-0.14 

0 

21+0.5 

■^-0.3 

6 Q+0-2 
-0.2 

-0.04 

296.77 

~0 

197.33 

~0 

99.44 

~0 

NLF2 

Q 1+0.13 

0 

00+5.6 

-3.2 

— 

-0.05 

39.2 

0.026 

21.45 

0.029 

17.75 

0.038 


Table 1. The best-fit parameters for an evolving LF (ELF) and non-evolving LF (NLF) for RGRBl and RGRB2 (with number 1 and 2 represent the two 
models), respectively. The statistical quality parameters of the models are shown in the last six columns. 


the table. For each model, we show the total and the goodness- 
of-fit Q (the probability to find a new;^^ exceeding the current one). 

In figure 3, we plot the observed and model-predicted PFD 
and RD. The left and right figures are the fitting results for RGRB1 
and RGRB2, respectively. The solid lines are for models with the 
evolving LFs, while the dashed lines are for models without evolu¬ 
tion. From this figure, we conclude our results as follows: 

1. For both GRB rate models, the evolving LFs (ELF) fit well 
with both PFD and RD, while the non-evolving LFs (NLF) fit the 
data poorly, which can be inferred from the values of the x^ and Q 
in Table 1. For RGRBl, the non-evolving LF (NLFl) model shows 
an excess of GRBs at high P and low redshift (left panels in figure 
3). For RGRB2 (NLF2), it also shows the same tendency as RGRB 1 
(right panels in figure 3). Therefore, we suggest that the LF should 
evolve with redshift and rule out the non-evolving LF for both GRB 
rate models. 

2. A non-negligible fraction of high-redshift GRBs may exist 
in the current Swift GRB sample. In our reduced redshift sample, 
the highest redshift is z = 6.32 (GRB 140515A). By using the best 
fit parameters provided in Table 1, we find nearly 23 GRBs with 
their redshifts higher than 6.5 for RGRBL Simulatively, nearly 24 
GRBs with redshifts higher than 6.5 for RGRB2. However, the truth 
is that we only observed 4 GRBs with possible redshifts higher than 
that (GRB 060116, z = 4 or 6.6; 080913, z = 6.44 or 6.7; 090423, 
z = 8.26; 090429B,z = 9.3). This result suggests that a large frac¬ 
tion of high-redshift GRBs (even with redshift higher than 10; e.g.. 
Tan et al. 2015) may exist in the BAT-triggered GRB sample. Un¬ 
fortunately, we can not measure their redshifts because of the in¬ 
strumental reasons. It is also possible that these high-redshift GRBs 
may origin from Pop III stars (e.g., Meszaros & Rees 2010; Toma 
etal. 2011). 

3. Except for the redshfit evolution of LF, an additional evolu¬ 
tion may exist between the GRB rate and SFR. If the LF and GRB 
rate model of RGRB2 represent the intrinsic ones, the additional 
redshift evolution is obvious while comparing the values of t be¬ 
tween the two models (r = 3.5);[52 for RGRBl and t = 0.8(;{5gg 
for RGRB2). However, it is hard to describe this evolution form 
as it may originate from many reasons, e.g., the evolution of cos¬ 
mic metallicity (Wang & Dai 2009; Li 2008), initial mass function, 
different properties of progenitor, their host galaxies, or their com¬ 
bination. The collapsar model explains the formation of GRB via 
the collapse of a rapidly rotation massive iron core into a black 
hole (Woosley 1993). However, there are two basic problems of 
single-star model in producing collapsars. First, the spin-down of 
the stellar core due to core-envelope coupling. Second, the Wolf- 
Rayet winds can slow the rotation. Fortunately, the rapid rotation of 
low-metallicity stars can overcome both problems (Langer & Nor¬ 
man 2006). Such speculation is confirmed by observation (Fynbo 
et al. 2003; Modjaz et al. 2008; Graham & Fruchter 2013; Trent! et 
al. 2013; Wang & Dai 2014). So we believe that it is easier for 



Figure 4. Redshift detection efficiency of BAT-triggered GRBs of Swift, 
with a total number of 244 GRBs with redshifts measured from afterglows 
and/or host galaxies. Top: Redshift detection efficiency for RGRBL The 
solid squares correspond to the detection efficiency of the total GRBs. 
The up-triangles (down-triangles) with dotted (dashed) line correspond to 
the detection efficiency of GRBs with peak fluxes > 10“^ erg s“* cm“^ 
(< 10“^ ergs^' cm“^). The solid fines are our fitting results. Bottom: The 
same as that of the top panel, but for RGRB 2. 

low-metallicity stars to produce GRBs than the stars with high- 
metallicity. If we incorporate this effect into the LF (as we have 
done in the model of RGRBl, where only the redshift evolution 
of LF is considered, but ignore the evolution of the GRB progen¬ 
itor properties), it will lead to a stronger redshift evolution of the 
LF, as the cosmic metallicity decreases with redshift. Futhermore, 
much more massive stars exist at high redshift, which will lead to 
more energetic GRBs and also increase the LF evolution. There¬ 
fore, the intrinsic redshift evolution of the LF should in the range 
of 0.8/;QQg < T < 3.5/;552 for reality. 

4.1 Redshift selection effect 

There are two major elements to affect the GRB redshift detection 
efficiency: one is the instrumental bias and the other is the redshift- 
dependent selection effect. For the first point, the localization and 
sensitivity of the follow-up afterglow observations have great im¬ 
pact on the detection probability of GRB redshifts. The quick re¬ 
sponse and high sensitivity making Swift detect the largest number 
and the highest redshift of GRBs. Assuming GRB optical after¬ 
glows decrease like power laws with exponent y -1, Swift will 
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detect GRB afterglows with 2.9 and 4.0 mag brighter than HETE2 
and BeppoSAX (Fiore et al. 2007), respectively. For the second 
point, the host galaxy dust obscuration plays a most important role 
in the measurement of GRB redshifts. It is suggested that a fraction 
of 30-35 percent of GRBs miss their redshifts for dust obscuration 
(Coward et al. 2013). In addition, the so-called “redshift desert” ef¬ 
fect can also lead to a fraction of nearly 20% of GRBs missing their 
redshifts in the range of 1.1 < z < 2.1 (Steidel et al. 2005; Fiore et 
al. 2007; Coward et al. 2013). The Malmquist bias arises because of 
the sensitivity of telescopes and instruments, which leads to a de¬ 
creasing redshift detection efficiency to higher redshift. Therefore, 
we only have 30% of GRBs with measured redshifts. 

In Figure 4, we compare the number distribution of the full 
redshift sample with our model predicted one in each redshift bin. 
The ratio could be considered as the redshift detection efficiency 
(solid squares). The error bars along the y-axis are the statistical er¬ 
rors (i.e., square root of the number in each bin AN - VN), which 
correspond to the 68% Poisson confidence intervals for the binned 
events (Gehrels 1986). The x-axis error bars simply represent the 
bin size. For both GRB rate models, the detection efficiency have 
the tendency of decreasing with redshift. Particularly, nearly 20% 
of GRBs miss their redshifts in the range of 1 < z < 2.5, which 
could be the evidence of the “redshift desert” effect. Furthermore, 
we divide the full redshift sample into two groups: one with peak 
fluxes P > 10“^ erg/ s/ cm^ (up-triangles with dotted line), and the 
other with P < 10“^ erg/ s/ cm^ (down-triangles with dashed line). 
It is obvious that the “redshift desert” effect mainly results from 
the low peak flux GRBs with P < 10“’erg/s/cm^. For GRBs 
with P > 10“’ erg/ s/ cm’, the detection efficiency just decreases 
with redshift and has no relation with the “redshift desert” effect. 
See from the figure, we could find more high (or less low) red¬ 
shift GRBs at low (high) peak fluxes, which may also indicate that 
high (low) redshift GRBs always have low (high) peak fluxes. For 
simplicity, we first fit the detection efficiency of RGRB1 with an 
exponential function multiple a Gauss function as 

rjAz) = 0.36 X x (1 - 0.41 x (12) 

Then we fit the sample with P > 10“’ erg/ s/ cm’ by multiplying 
an constant with equation (12) but without the Gauss fraction, and 
the fitting result could be expressed as 

j/hCz) = 0.16xe““=^’**‘’. (13) 

For the detection efficiency of P < 10“’ erg/ s/ cm’, we obtain it by 
subtracting equation (13) from (12), which is r]\{z) = pAz) ~ flh(z)- 
It seems that p](z) also fits the reuslt quite well, which could also 
support our result that the “redshift desert” effect mainly comes 
from low flux GRBs. For RGRB2 model, we go through the same 
steps and the fitting results are as follows, 

p[(z) = 0.33 X X (1 - 0.3 x g-fe-i T’/oo’) (U) 

and 

^/[.(z) = 0.15xe°’“^’‘’«. (15) 

Finally, we obtain p[(z) = p[(z) - t 7 f,(z), which also fits the obser¬ 
vation quit well. The fitting results are shown as the solid lines in 
Figure 4. The goodness of the fits at the highest redshift are greatly 
affected by the GRB number. Here we should mention that the “red¬ 
shift desert” effect exists in both GRB rate models, particularly in 
the dimmer GRB sample. However, if this is caused by the statis¬ 
tical fluctuations, the “redshift desert” effect will disappear. Then, 
the redshift detection efficiency will be nearly a constant at z < 2.5 
and decrease at higher redshift. The problem is that the fluctuations 


only exist in the dimmer GRBs, but not in both dimmer and brighter 
GRB samples. The so-called “redshift desert” is a redshift region 
(l<z<2.1) where it is difficult to measure absorption and emis¬ 
sion spectra of the sources. As redshift increases beyond z ^ 1, the 
strong emission lines from the host galaxies (e.g. the [O III]/t/t4959, 
5007, [O II]/I3727 line, Ha, H/l) are shifted outside of the typi¬ 
cal interval covered by optical spectrometers (x 3800 - 8000A) at 
z > 1, while Lyman-a enters the range at z ~ 2.1. This becomes 
more serious for faint sources, as the spectra of the faint sources 
always have small signal to noise ratio. Therefore, we believe that 
dimmer GRBs are more responsible for the “redshift desert” effect, 
which can not be caused by the statistical fluctuations. 


5 SUMMARY AND DISCUSSION 

In this work, we utilized both the peak flux and redshift distribu¬ 
tions to constrain the LF. We don’t use the pseudo redshifts in¬ 
ferred from the empirical luminosity relations, because these rela¬ 
tions could be greatly affected by the observational bias, e.g., trig¬ 
ger efficiency (Shahmoradi & Nemiroff 2011). By jointly fitting 
the PFD and RD, we found that the non-evolving LF fits the data 
poorly, while the evolving LF fits data well. The evolution of the 
break luminosity can be expressed as Lt oc (1 -t-z)’^, with t = 3.5/q2 
and T = 0.8/QQg for RGRBl and RGRB2, respectively. The strong 
and weak redshift evolution of the LFs for two GRB rate models 
might imply an additional redshift evolution between GRB rate and 
SFR. The strong redshift evolution of LF for RGRB 1 is based on 
the assumption that only LF evolves with redshift, which is not real 
if GRB progenitors evolve with redshift. If we take these effects 
into account, the redshift evolution of the LF may become weaker 
for RGRBl (e.g., t < 3.5/[j2). However, it is impossible to distin¬ 
guish this additional evolution under the current situation, unless 
the number of GRBs with measured redshifts is large enough to 
satisfy the LF fittings in different redshift ranges (Tan et al. 2013). 

The local GRB rates are R(0)i = 0.86/[5gg yr“* Gpc“’ and 
R(0)2 = 0.54/“;’^ yr“' Gpc“’ for RGRBl and RGRB2, respectively, 
which are consistent with the previous works (e.g., Schmidt 1999, 
2001; Guetta et al. 2004, 2005; Liang et al. 2007; Wanderman & 
Piran 2010; Gao et al. 2011; Lien et al. 2014). Furthermore, we 
suggest that a large fraction of high-redshift GRBs (z > 6.5) may 
have been detected by Swift satellite but without measured red¬ 
shifts. These high-redshift GRBs may be produced by Pop III stars 
(e.g., Meszaros & Rees 2010; Toma et al. 2011; Tan et al. 2015). 

To reduce the large biases introduced by the complex trigger¬ 
ing criterion of Swift, we use the recent peak flux efficiency func¬ 
tion for the triggering criterion instead of the constant photon flux 
(Howell et al. 2014). The constant photon flux triggering criterion 
could be problematic while estimating the contribution of the dim 
and bright GRBs as shown in figure 2. The introduction of a depen¬ 
dence of L on Ep has little effect on our results, and the Yonetoku 
relation is considered to be tight enough. 

To reduce the redshift selection effect, only GRBs with red¬ 
shifts measured from afterglows or from both afterglows and host 
galaxies are considered. After deriving the best-fit parameters, we 
compare the observed full redshift distribution with our model 
predicted one. Our results show that the redshift detection effi¬ 
ciency decreases with redshift slowly in general, and nearly 20% 
of GRBs missed their redshifts in the “redshift desert” range of 
1 < z < 2.5. The reason for the slow decreasing tendency of the de¬ 
tection efficiency might be the combined effects of the dust extinc¬ 
tion and Malmquist bias, as the higher (lower) redshift GRBs al- 
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ways have less (more) dust contents but lower (higher) peak fluxes 
(e.g. Bouwens et al. 2009; Perley et al. 2009; Zafar et al. 2011; 
Rossi et al. 2012). By dividing the GRB sample into the dim and 
bright groups, we find that the “redshift desert” elfect is mainly 
introduced by the low-flux GRBs, which is a useful result for the 
completeness of the sample selection in future works. 
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